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DEGRADATION AND HEALING IN A GENERALIZED NEO-HOOKEAN 
SOLID DUE TO INFUSION OF A FLUID 

SATISH KARRA AND K. R. RAJAGOPAL 

Abstract. The mechanical response and load bearing capacity of high performance poly- 
mer composites changes due to diffusion of a fluid, temperature, oxidation or the extent 
CN , of the deformation. Hence, there is a need to study the response of bodies under such 

degradation mechanisms. In this paper, we study the effect of degradation and healing due 
to the diffusion of a fluid on the response of a solid which prior to the diffusion can be 
described by the generalized neo-Hookean model. We show that a generalized neo-Hookcan 
solid - which behaves like an elastic body (i.e., it does not produce entropy) within a purely 
mechanical context - creeps and stress relaxes when infused with a fluid and behaves like 
a body whose material properties are time dependent. We specifically investigate the tor- 
sion of a generalized neo-Hookean circular cylindrical annulus infused with a fluid. The 
equations of equilibrium for a generalized neo-Hookcan solid are solved together with the 
convection-diffusion equation for the fluid concentration. Different boundary conditions for 
the fluid concentration are also considered. We also solve the problem for the case when the 
q diffusivity of the fluid depends on the deformation of the generalized neo-Hookean solid. 
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>v 1. Introduction 

^ . 

Q_i! As elastic bodies are incapable of producing any entropy, this would be the proper defini- 

tion of an elastic body within a thermodynamic context. Elastic bodies also cannot stress 
relax, i.e., when the strain is kept constant in time the stresses within the body cannot change 
with time. However, the stresses in a body that is initially elastic, which subsequently un- 

CH ! dergoes chemical reactions due to interactions with the environment or which is subject to 

the effects of electro-magnetic radiation such as ultra-violet rays, could change with time 
and the body's response ceases to be that of an elastic body. The chemical reactions or the 
interactions with the environment invariably produce entropy. However, it could happen 
that when the chemical reactions cease and the body is isolated from the environment the 
body ceases to produce entropy, that is, it becomes a different elastic body. The body, due 
to its exposure to the environment can undergo deterioration or enhancement with respect 
to its load carrying capacity or other properties. While moisture diffusion in a polymer can 
cause degradation of the body in that its load carrying capacity goes down, a body such 
as biological matter can be strengthened due to a drug that is being injected. Though the 
stress in the body might decrease when the strain is held constant in bodies that are under- 
going degradation, this phenomenon is quite different from the stress relaxation observed in 
viscoelastic solids (Christensen [5], Wineman and Rajagopal |32j). 

There is need for a good understanding of the degradation of materials as this is relevant 
to a plethora of applications. For instance, the mechanical properties of high performance 
polymer composites (like polyimides) used in hypersonic vehicles vary due to the effect of 
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high temperature, diffusion of moisture and the subsequent oxidation. Chen and Tyler [1] 
have recently shown that polymers can also degrade due to deformation. In this short paper 
we are interested in studying the response of a body whose properties are changing due to 
the presence of a fluid, the extent of the change depending on the concentration of the fluid. 

There has been considerable interest in the damage and degradation of polymeric solids and 
polymer composites. However, most of these studies appeal to ad hoc constitutive equations. 
Also, many studies appeal to the notion of a hygrothermal expansion coefficient (see Bouadi 
and Sun [I], Cai and Weitsman [3]) and such an approach would not be appropriate when 
large deformations due to swelling are involved. A general thermodynamic framework has 
been developed to describe damage in composites by Weitsman [29]. This work introduces 
the notion of a damage tensor which is essentially a tensor internal variable. While the 
study is very detailed and has a thermodynamic basis, the theory however involves several 
material moduli that depend on as many as 32 invariants in the case of transversely isotropic 
materials, making it impractical to put the theory to use as it is impossible to develop an 
experimental program to determine the numerous material functions that characterize the 
body. Moreover, such theories lead to intractable mathematical equations. While there have 
been several other studies concerning the diffusion of moisture in composites (e.g., Snead 
and Palazotto [23] , Lee and Yen [13] , Bouadi and Sun [2] , and Kardomateas and Chung [II] ) , 
they are primarily parametric studies in which the material moduli are allowed to depend 
on the moisture content according to some pre-assigned manner and not as a consequence of 
a reasonable convection-diffusion equation; that is the coupled problem for the deformation 
of the composite and the diffusion of moisture is not solved simultaneously. In general, the 
diffusivity depends on the temperature, moisture content, strain, and stress. In fact, the 
diffusivity can depend on the history of these quantities. Weitsman [3U] and Roy et al. [2T| 
have also addressed damage due to diffusion but not in the manner advocated here. 

Using ideas in multi-network theories for polymeric materials (see Tobolsky and Andrews 
[26] . Tobolsky et al. [27]), Wineman and Rajagopal [31] and Rajagopal and Wineman [T9] 
developed a theory applicable to the large deformation of polymeric solids that exhibit scisson 
and healing and permanent set, within a purely mechanical context. This work was extended 
in a series of papers by Huntley et al. [TU], [5], [Sj to deal with a variety of deformations 
involving elastomeric solids that undergo permanent set. The scission and healing that 
takes place can be viewed within the context of deformation induced damage and due to 
cross-linkings that take place. These works within a purely mechanical context were later 
generalized to include thermal effects (see Wineman and Shaw [33], Shaw et al. [23]). 

A general thermodynamic framework which takes into account chemical reactions (with 
chemical kinetics), diffusion and thermal effects needs to be put into place in order to develop 
constitutive relations to study the degradation of materials (like polyimides) due to chemical 
reactions. As a first step towards such a goal, in this paper, we shall first solve the problem 
of torsion of a degrading body, which when there is no degradation taking place responds like 
a generalized neo-Hookean body. We also look at healing (or strengthening) of a generalized 
neo-Hookean body when a fluid is infused. We assume the body to be a cylindrical annulus 
of finite length. We shall study the torsion of a cylindrical annulus through which a fluid is 
infusing. We will assume that the infusion of the fluid is radial and thus the degradation or 
healing takes place radially. We introduce a parameter that is a measure of the degradation 
or healing which in virtue of the diffusion being radial also varies radially. The material 



parameters are assumed to be functions of the variables that quantify the degradation or 
healing. We need to then solve the convection-diffusion equation that governs the diffusion 
of the fluid in tandem with the equilibrium equations for the torsion problem; we look at how 
the moment varies with time when the angle of twist of the cylinder is kept constant in time, 
and how the angle of twist of the cylinder varies when the moment applied to the cylinder is 
kept constant. We find that the moment that is needed to maintain the angular displacement 
decreases with time, that is the body stress relaxes, when the material is degrading, but as we 
observed earlier such a decrease of stress is very different from the stress relaxation observed 
in viscoelastic bodies. The stresses can decrease in the body due to a variety of reasons 
such as degradation of properties of the body due to chemical reactions, aging, etc., and it 
is important to recognize the reason for the "stress relaxation". Rajagopal and Wineman 
[20] have studied stress relaxation due to aging and they find that in marked contrast to the 
stress relaxation observed in viscoelastic bodies, the decrease in stress is dependent on the 
geometry of the body. This aspect related to stress relaxation, namely its dependence on the 
geometry is what sets stress relaxation in viscoelastic materials apart from stress relaxation 
that manifests itself due to most degradation theories. In order to highlight this difference, 
Rajagopal and Wineman [20] considered the torsion of a viscoelastic cylinder that ages. 
They found that the stress relaxation can be split into two parts, one that is a consequence 
of the body being viscoelastic, which is independent of the size of the specimen, and another 
part due to the aging of the cylindrical specimen, this being dependent on the radius of the 
cylinder undergoing torsion. We also find that the angular displacement undergone by the 
body increases with time for the applied moment, when the material is degrading. When one 
looks at healing due to diffusion of a fluid, we find that that the moment needed to maintain 
a given angular displacement increases with time, whereas the body's angular displacement 
decreases with time when one applies a constant moment. 

We will now turn to a discussion of the response characteristics of the undamaged elastic 
solid, namely the power-law neo-Hookean solid. The generalized (or power-law) neo-Hookean 
elastic model (see Knowles [12]) allows for softening and stiffening under simple shear when 
the power-law parameter (n) is lesser or greater than unity, respectively. Softening and stiff- 
ening which are seen in real materials cannot be explained using the classical neo-Hookean 
model while some of it possibly can be explained in terms of the generalized neo-Hookean 
solid. In any event, one should recognize that such models are merely caricatures of reality 
and in this study we are mainly interested in obtaining some understanding of the degrada- 
tion due to the infusion of fluid. On setting the power-law parameter (n) to 1, one obtains 
the classical neo-Hookean model. Knowles [12] studied the anti-plane shear of a power-law 
neo-Hookean body that has a crack, and noted that for n > |, the equation of equilibrium for 
anti-plane shear is always elliptic and that the ellipticity is lost when n < ~. Hou and Zhang 
[7] have studied the stability of a power-law neo-Hookean cylinder under axial stretching 
and constant radial traction. The power-law model has also been used in several subsequent 
studies (see Saccomandi [22]). 

Rajagopal and Tao [TS] studied inhomogeneous deformations in a wedge of a generalized 
neo-Hookean material numerically, and later, Mcleod and Rajagopal [H] re-investigated 
the same problem with a view towards establishing existence of solutions to the governing 
equations. Tao and Rajagopal [25] have analyzed the problem of inhomogeneous circular 
shear combined with torsion for a genearalized neo-Hookean material and have obtained 



exact solutions for certain values of the material constants. This study was followed by 
an analysis by Zhang and Rajagopal [34] concerning steady and unsteady inhomogeneous 
shear of a slab and a cylindrical annulus. The reason for citing these studies is that in 
all of [TSJ dU [221 EI], "boundary layer" type solutions (that is, close to the boundary the 
deformations are inhomogeneous while the deformations are mainly homogenous in the inner 
core region of the body) have been observed for certain values of the material parameters 
and also the strain is concentrated in these boundary layers in that the strain gradients are 
very large in these layers. Thus, one would expect damage, degradation and failure to occur 
in these boundary layers. A systematic method to obtain approximate equations within the 
boundary layer of a generalized neo-Hookean solid similar to the boundary layer theory in 
fluid mechanics has been discussed by Rajagopal [17]. To show the efficacy of this method, 
the approximate boundary layer equations were solved for the circumferential shear problem 
and were compared to the full solution. 

While there have been several studies concerning generalized neo-Hookean materials within 
the context of a purely mechanical setting, as noted above, there is no study concerning the 
degradation or enhancement of the properties of generalized neo-Hookean materials, which 
would be important as materials such as elastomers for which the generalized neo-Hookean 
model is used, degrade or heal due to the diffusion of moisture or other chemicals. 

Recently, there have been a few investigations of the deformation of an elastic solid that 
is undergoing degradation due to the influence of a diffusant. Muliana et al. [15] analyzed 
the response of a composite cylinder that is undergoing degradation due to the diffusion of 
a fluid, and Darbha and Rajagopal [H] studied unsteady motions of a slab through which a 
fluid is diffusing. They investigated the unsteady shear of an infinite slab of finite thickness 
and a cylindrical annulus of infinite length under degradation due to diffusion of a chemical 
species. In both the studies the body was assumed to be a linearized elastic body and the 
material moduli were assumed to vary linearly with the concentration of the diffusing species. 

Rajagopal [16] assumed the shear modulus of a generalized neo-Hookean material to de- 
pend on temperature and solved the problem of inhomogeneous shear of an infinite slab. 
Boundary layer type solutions were obtained based on the nature of the boundary con- 
ditions for the temperature. Although, this paper was not written within the context of 
understanding degradation due to temperature, one can extend this work by assuming that 
all the material moduli depend on temperature in a certain fashion and study the effect of 
temperature on the response of the body. It is expected that boundary layers would develop 
due to degradation and the failure of the body would be determined by the stresses in this 
boundary layer. 

The current paper is organized as follows. In section [21 the problem of torsion of a finite 
cylindrical annulus comprised of a generalized neo-Hookean material through which a fluid 
is diffusing, is set up. In section El we specifically assume that the degradation or healing is 
due to diffusion and also assume that material moduli vary linearly with the concentration of 
the diffusant, and obtain the solutions to the convection-diffusion equation for the diffusing 
species as well as the equilibrium equations. We also obtain an expression for the moment 
that is applied in terms of the the angular displacement and the material parameters. In 
section [H the results are discussed in detail, followed by final conclusions in section [51 



2. Torsion of a Cylindrical Annulus undergoing Degradation 



Let kr and n t denote the reference and the current configuration of a body, respectively. 
The motion Xk r is a one-one mapping that assigns to each point X e Kr, a point x e n t , at 
a time t, i.e., 

x = x Kfl (X,t). (1) 

The gradient of motion (or the deformation gradient) F is defined by 

dXn R 



with the velocity v defined as 



F:= 



v :- 



<9X ' 

9 Xk r 

dt ' 



(2) 



(3) 



Let (R, 0, Z) and (r, 6, z) be cylindrical co-ordinates in the reference and current configura- 
tion, respectively. Consider a cylindrical annulus of height H with inner radius Ri and outer 
radius R under torsion, whose motion is given by 

r = R, 9 = Q + f(z t t), z = Z. (4) 

The deformation gradient F and the left Cauchy-Green stretch tensor B associated with (J3J) 
are given by 



1 
1 rf z 
1 



B := FF J 



1 

\ + {rf z f rf z 
rf z 1 



(5) 



df(z t) 
Thus, the first invariant of B, I± = tr(B) = 3 + (rf z ) , where f z := — — - — , and tr(.) is the 

trace of a second order tensor. 

The Cauchy stress tensor T in an incompressible generalized neo-Hookean material, is 
given by 

n-l 



-pi + (J, 



1 + - (h - 3) 

n 



B. 



(6) 



where —pi is the spherical part due to the constraint of incompressibility, and fj, [1 + - (I± — 3 
is the generalized shear modulus, /i being the shear modulus at zero stretch. In general, the 
degradation or healing of the generalized neo-Hookean material can be caused by diffusion, 
temperature, electromagnetic radiation, etc. We shall denote the variable that is a measure 
of the degradation or healing by a, i.e., when degradation due to diffusion is considered, 
a would be the concentration of the diffusing species; a would be temperature if we have 
degradation or healing due to temperature, etc., with an appropriate governing equation for 
the variable. We shall refer to a as the degradation/healing parameter. We shall assume 
that the material moduli are functions of the degradation/healing parameter i.e., n = n(a>), 
fi = /i(a), b = b(a), which leads to changes in the response characteristics of the body. In 
view of the geometry of the body, and the assumed form for the deformation field, we shall 
further assume that the degradation/healing parameter varies only with the radius and time, 
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n-l 



i.e., a = a(r,t). Then, the material moduli would be functions of the radius and time and 
(ED) reduces to the form 



b(r,t) ( f 
n(r,t) V Jz 



pI + $(r,z,t)B, (7) 

n(r,t)-l 



where $(r, z, t) := fi(r, t) 

We shall consider the diffusion to be reasonably slow and that the inertial effects in the 
solid can be neglected. It is important to bear in mind that the temporal effects are not 
being ignored. One could choose to view time as a parameter. Neglecting the body force, 
the balance of linear momentum reduces to the equations of equilibrium, 

div T = 0, (8) 

where div(.) denotes the divergence operator in the current configuration. We assume that 
the concentration of the diffusing species (c) is governed by the following convection-diffusion 
equation 

dc 

— + div (cv) = div (D grad(c)) , (9) 

L/C 

where D is the diffusivity and in general it could depend on the deformation as well as the 
concentration, grad(.) is the gradient based on the current configuration. In this study, we 
assume that the diffusion is very slow so that the velocity of the fluid that is diffusing can be 
neglected. The second term on the left side of (J9]) involves the derivative of the velocity and 
of course it is possible that even if the velocity is small its derivatives need not be small. We 
shall however assume that the derivatives of the velocity are small and thus we neglect the 
second term in ((9]). We shall assume that D depends on the deformation, thus the equation 
governing the diffusion of the fluid is coupled with the balance of linear momentum. Since 
the concentration is only a function of the radius r and time t, equation (Q reduces to 

dc(r,t) 1 d f Dr . 9c(r,t) \ (10) 



dt r dr \ dr 

Next, we shall document the balance of energy for the generalized neo-Hookean solid and 
fluid as follows (see Truesdell et al. [28]): 

dp 1 
tf = TW - divq 4 + pV + E\ i = solid, fluid, (11) 

dt 

where e\ q\ r % are the specific internal energy, heat flux, radiant heating associated with 
the i-th component and E l is the energy supplied to the z-th constituent from the other con- 
stituents. We shall ignore the energy equation associated with the fluid and the contributions 
due to the interactions between the fluid and the solid {E s ) in the energy balance for the 
solid. We shall however later on incorporate the energy associated with fluid indirectly by 
assuming that the internal energy of the solid depends on the concentration of the diffusing 
fluid. This is clearly a crude approximation. In such a case, the energy equation is merely, 

pe = T.L - divq + pr, (12) 

where e is the internal energy of the solid, L is the velocity gradient of the solid, q is the 
heat flux, and r is the specific radiant heating. Further assumption that the internal energy 
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of the solid depends on the temperature (T), concentration of the diffusing fluid and the 
deformation gradient of the solid, i.e., 

e = e(T,c,F), (13) 

leads to 

de ■ de de ^ T t m t , N 

p df T + P ^ C + P 9F F L = T L ~~ dlVq + Pr ' (M) 

where (.) T is the transpose of a second-order tensor. We shall define the specific Helmholtz 
potential as \P — e — Trj, where rj is the specific entropy. Then ( Tl4l) reduces to 

p %f + „*V.L - P T fl-F T .L + P ^c - pT^-c = T.L - divq + pr. (15) 

H dT ' dF ' dTdF ' dc ' dTdc H H K J 

dty d^ff 

where we have also used that 77 = — — — . Now, if we set T = — 71T + p-^F T , where ir is the 

Lagrange multiplier due to the constraint of incompressibility (given by tr(L) = 0), then in 
the absence of radiation along with the assumption of Fourier's relation for heat conduction, 

q = —k grad(T), (16) 

( 1T5|) reduces to 

d 2 ff - &$ d 2 ff 

pCvT - pT^^F^L + p—c - pT— c 6 = div (k grad(T)) , (17) 

de d 2 if> 

where C v = -7-= = —T——- is the specific heat capacity and k is the heat conductivity, 

oT 0T 1 

both of which could depend on the deformation, temperature as well as concentration of the 

diffusing fluid (in general, it could depend on the degradation/healing parameter a). We 

also note that for a generalized neo-Hookean body, the specific Helmholtz potential is given 

by 



-&{K<H"- 1 } 



In general, along with infusion of a fluid, the body could be subject to high temperature 
which could cause additional degradation, and so the material parameters could depend on 
temperature as well i.e., n = n(c,T), p = p(c,T), b = b(c,T). In such a case, one ought to 
solve flSJ), ( ITU|) . (TT7|) . simultaneously along with (ITSjl . In what follows, we shall ignore the 
thermal effects and study the problem of degradation or healing only due to the diffusion of 
a fluid. 

3. Degradation and Healing due to Diffusion 
It follows from flZJ and (JHJ that 

— (-p + $(r, z, t)) - r^{r, z, t) {f z f = 0, (19) 

-qq(-P + Hr, z, t)) + — ($(r, z, t)rf z ) = 0, (20) 

^(-p + $(r,z,t)) = 0. (21) 



In solving ( TT9|) . ( 120]) and (l2~il . time £ is treated as a parameter i.e., we are solving the balance 
of linear momentum flH]) at every instant of time assuming that the problem is quasi-static. 
Next, from (ED), it follows that 



-p + $(r,z,t)=g(e,r). (22) 

Taking the derivative of (fT9|) with respect to 9, and using (122]) . we obtain that 

8 fdg(r,i 



8r \ 89 
dg{r, 9) 
89 



0. 



h(9), 



and therefore (1221 . (123]) reduce to 

90 



(-p + $(r,z,t)) = h(9), 



where h is a function of 9. From ( 120]) and (T24|) . we obtain that 

-^ (-P + $(r, 2, *)) = r— ($(r, z, t)r/ z 



£. 



(23) 



(24) 



(25) 



9p 



where -E is a constant. On the basis of periodicity, we can assume that — = 0, and hence 

89 
from ( 125]) it follows that E must be zero. Using the definition of $ in (125]) . we obtain that 



where f z 



1 + - (rf z ) 2 
n 



d 2 f(z,t) 



-in-l 



nn-D (rh? 



n 



b o 

1 + - (rf z ) 2 
n 



n-2 



fzz 



E 
fir 2 



0, (26) 



8z 



. Thus, (126]) reduces to either 
fz, = 0, 



or 



1 + - {rf z 

n 



Equation (T27]) implies that 



and (1281) reduces to 



n-l 



26( - 1) (r /,)' 



/? 



1 + - (r/,) 2 
n 



-i n-2 



f(z,t) = C 1 (t)z + C 2 (t) 
b(2n - 1 



1 



n 



{rf z 



0. 



(27) 



(28) 



(29) 



(30) 



which has no solution. Hence, ( 129]) is the solution to (jHJ). Interestingly, the solution ( 129]) 
does not depend on the material parameters n, b or \x. 

If we assume that the angle of twist at z = is zero, i.e., if we assume that the bottom of 
the cylinder is fixed, then we obtain that 



f{z,t)=4>{t)z, 



(31) 



where ip{t) is the angle of twist per unit length of the cylinder which is the time dependent 
generalization of the classical torsion solution. The twisting moment is given by 



M(t) 


rRo 
= 2n T z gr 2 dr 
Jr 1 






rRo 

= 2ir ti(c(r,t)) 

JR 1 


I ! K<r,t)) 
n(c(r,t)) 



- n(c(r,t))-l 

{rip(t)) 2 ip(t)r 3 dr. (32) 



In our work, we shall enforce the following two sets of initial and boundary conditions for 
the concentration of the diffusing species: 
Case 1: 



c(r, 0) = 0, R % <r< R , (33) 

3c 

— {Ri,t)=0, Vt>0, (34) 

c(R ,t) = 1, Vt >0. (35) 

That is, initially, the body is assumed to be in its virgin state and there is no diffusant in the 
body. Also, the boundary condition (134)) implies that the gradient of the concentration is zero 
at the inner boundary of the annular cylinder. This can be achieved by having some sort of 
a membrane on the inside of the annulus which is impermeable to the diffusing species. The 
boundary condition (I3"5"|) means that the outer cylinder is always exposed to the diffusant. 
Case 2: 

c{r, 0) = 0, Ri<r < Rq, (36) 

c(Ri,t) = co, Vt>0, (37) 

c(R ,t) = 1, Vt>0. (38) 

Here, as opposed to case 1, the inner boundary of the annular cylinder is held at a con- 
stant diffusant concentration (as reflected by f )37|) ). By constructing a mechanism which 
continuously removes the diffusing species from the inside of the cylindrical annulus, a con- 
stant concentration boundary condition can be maintained. For example, if the diffusing 
species is moisture, then by blowing air inside the annulus continuously, one can control the 
concentration of the moisture. For convenience, we shall set cq to zero for this case. 

In addition, we shall assume that the material moduli change in a linear fashion due to 
the diffusion, that is the material parameters change with c as follows: 

H = fJ, ± fJtic, < fjt! < Ho, (39) 

b = b ± he, < 6i < b , (40) 

n = n ± riic, < n\ < n , (41) 

with the minus sign chosen when degradation is considered and plus sign for healing. The 
restrictions on the ranges of /ii, b±, n\ are being enforced to ensure that /x, b, n are positive, 
in the case of degradation. 
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Before proceeding with the solutions, we render the governing equations dimensionless. 
The non-dimensionalization is carried out in the following manner. We define 



t 



t r z - 

— ,r = — ,z = — ,ip 



t 



Rn 



W 



6q /i 



(42) 



where to, #o are characterstic time and angle. Then (132]) and (fTOj) reduce to 



MH 



27tR%fi , 

M 



l> 



-1 ra-1 



n H z 



^3^r- 



ipr dr 



(43) 



dc(r,i) 1 d 
di f df 




(44) 



where /i = 1 ± \i\c, with \x\ 



Mi 

Mo' 



The initial and boundary conditions for c after non-dimensionalization reduce to 
Case 1: 

c{f, 0) = 0, fj < f < 1, 

dc - 

— (f h t) = 0, Vt>0, 

c{l,t) = 1, W>0. 



Case 2: 



, _ R% 

where r, = — . 



c(r, 0) = 0, r, < r < 1, 
c(f i; f) = 0, Vt > 0, 
c{l,t) = 1, Vt >0, 



4. Discussion of Results 



(45) 

(46) 
(47) 

(48) 
(49) 
(50) 



We now proceed to solve the non-dimensionalized problem, numerically. Figure (I lap shows 
the solution to the convection-diffusion equation ( 1441) for D = 0.01, which has been solved 
using pdepe solver in MATLAB. We notice that after a certain time, the concentration profile 
reaches a steady state which implies that no further healing or degradation of the body takes 
place. For a given value of ip, the values of concentration thus obtained at various radii and 
times were used to numerically integrate (143|) using composite trapezoidal rule to find M. 
On the other hand for a given value of M, ip was calculated using a combination of the 
bisection method and the composite trapezoidal rule of intergration on (143|) . The numerical 
values chosen for the non-dimensional quantities are shown in the figures. 

The solution depicted in figure (fib"]) was obtained by setting n to 1 along with n\ to 
when degradation is assumed. Under these conditions the model reduces to a neo-Hookean 
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model as n = 1; furthermore, note that p,\ = corresponds to case when the neo-Hookean 
body is neither degrading nor healing. Notice that as Ji\ increases, the stress required to 
maintain the angular displacement decreases, i.e., the higher the degradation of the material, 
the greater is the stress relaxation in that less moment is necessary to maintain the angular 
displacement. Of course, since the moment is related to the shear stress through (132 p . we can 
call this stress relaxation due to degradation. In figures (CQ), (J2J) the non-dimensional moment 
is portrayed as a function of non-dimensional time by varying /i 1; b\ and n\ along with other 
quantities as indicated. In all these cases, it is seen that as the degradation increases, less 
moment is needed to maintain the deformation, as is to be expected. Next, we assume that 
the material is healing as the fluid diffuses (see figure (J3J)), and the non-dimensional moment 
is observed as a function of time for a fixed value of non-dimensional angular displacement. 
As the body heals with time, the moment required to maintain the deformation increases. 
Furthermore, as seen from figures (CQ), (|2j), (jSJ) the non-dimensional moment reaches a steady 
value to maintain the angular displacement since the concentration of the diffusant reaches 
a steady value after which there is no further degradation or healing. 

Figure 01]) shows that the non-dimensional angular displacement of the cylinder increases 
with time for an applied non-dimensional moment when the generalized neo-Hookean body is 
degrading i.e., the body "creeps" for sometime, and then the angular displacement reaches a 
steady value. This is very different from the creep in a viscoelastic solid, wherein the angular 
displacement continuously increasing with time when an external moment is applied. This 
is one way of determining if the creep undergone by a body is either due to the viscoelastic 
nature of the body or due to degradation. Now, when one considers healing, the non- 
dimensional angular displacement of the body decreases with time as shown in figure (jSJ) for 
an applied moment, and then remains steady. 

As mentioned in the introduction, Rajagopal and Wineman [20] showed that the stress re- 
laxation due to aging depends on the material geometry, and this characteristic differentiates 
the stress relaxation due to aging/ degradation from the stress relaxation due to viscoelas- 
ticity. To illustrate this phenomenon in our work, results were obtained at the ratio of the 
inner radius of the annulus to the outer radius (fi) being 0.35 and 0.75 (see figure ©), and 
it is seen that not only stress relaxation but also creep due to degradation depends on the 
geometry. Comparing figures (jBaj) and (J6bl . figures (|5cj) and (J6dj) . one can also see that the 
steady state values are attained at a faster rate when fj = 0.75. This is because the annulus 
has a smaller thickness for this case and so the concentration reaches steady state faster. 

To see how the stress relaxation and creep depend on the diffusivity, we changed the 
non-dimesional diffusivity (/)) from 0.01 to 0.1. As can be seen in figure ([7]) the moment 
relaxes faster as the diffusivity increases, and the angular diplacement undergone by the body 
increases faster as the diffusivity increases. This is because higher diffusivity means that the 
concentration at any given location increases faster as it tends towards the steady state value, 
and as the material parameters decrease or increase (based on degradation or healing) in 
value with increasing concentration, the material degrades or heals faster. Also, notice that 
the steady state non-dimensional moment and the non-dimensional angular displacement 
values are the same for the three values of the non-dimensional diffusivity. This is due to 
the fact that the steady state solution for concentration of the diffusant is the same in all 
the cases, which is c(r,t) — 1, for t sufficiently large so that steady state is reached. 



n 



Next, we shall look at how the results vary when the diffusivity depends on the strain; we 
choose the Almansi-Hamel strain as the measure of the non-linear strain. This is to capture 
the fact that the pore structure in the body depends on the strain and hence will lead to 
a change in the diffusivity. However, apriori, it is not clear whether the diffusivity has to 
increase or decrease. In this study, we shall assume that the diffusivity increases with the 
strain. 

Now, the Almansi-Hamel strain is given by 



e=i(I-B"') 










—rip 



(51) 



The Frobenius norm of e is given by 
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e T e) 
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J_ q 2 ( r ^/) 4 + _ g ( r ^) 2 . 



(52) 



We shall assume that the diffusivity D depends on the Almansi-Hamel strain in the following 
manner 

D = Do + (Ax> - Do) (1 - e^ |e|1 ) . (53) 

Figures ( !8cl) . ( |8dl) compares the stress relaxation and creep for the case with constant 
diffusivity against the case when the diffusivity depends on the Almansi-Hamel strain with 
the initial and boundary conditions considered are given by (HoTHTl) . The stress relaxation 
and creep are faster when diffusivity is assumed to be a function of the Almansi-Hamel 
strain, as the diffusivity is higher in this case and hence the degradation is faster. 

Next, we shall consider the situation when we have the initial and boundary conditions 
given by (J48H50]) . Here too, the stress relaxation and creep are faster when the diffusivity 
depends on the Almansi-Hamel strain is assumed (see figures (I9cl) . (|9d|) ). However, the steady 
state values of the non-dimensional moment and non-dimensional angular displacement when 
the diffusivity is maintained constant is not the same as that when the diffusivity depends 
on the Almansi-Hamel strain. One can explain this by looking at the steady state solutions 
for the convection-diffusion equation in figure fllOl) . Unlike the situation corresponding to 
(H5H171) . the steady state solutions for the concentration are not the same and the steady 
state concentrations when the diffusivity depends on the strain are higher than the steady 
state values when the diffusivity is constant. This induces the stresses to relax to a lower 
value in fl9cl) and the angular diplacement to reach a higher value in (I9d|) . These values are 
also attained faster as the steady state solution is reached faster when the diffusivity depends 
on the strain. 



5. Conclusions 

Using the generalized neo-Hookean solid as a vehicle to illustrate the effects of degradation 
of a body that is initially elastic we have shown that degradation due to the infusion of a 
fluid can cause stress relaxation and creep until a steady state is reached. In contrast to a 
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viscoelastic body which creeps continuously upon application of a load, the body considered 
here stops to creep after a certain steady state value is attained. This is one important 
characteristic that can be used to differentiate the creep due to the viscoelastic nature of 
a body, and the creep due to degradation. We have also considered the case of healing or 
strengthening of a generalized neo-Hookean body due to the infusion of a fluid. One can 
easily extend this work to other forms of degradation in a generalized neo-Hookean material 
or for that matter in any elastic body whose properties change due to the infusion of a 
fluid, by providing an appropriate equation that governs the evolution of the degradation 
parameter. 
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Figure 1. (a) Solution to the convection diffusion equation, (b) Non- 
dimensional moment (M ) as a function of non-dimensional time it) for various 
values of Jx\ starting from to 0.5 in increments of 0.1 for the degradation case. 
Values chosen were Ri/R = 0.5, q — 1, if) — 1, D — 0.01, b = n = 1, 6 X = 0, 
ni = 0. This corresponds to the neo-Hookean model since n — 1. (c) Non- 
dimensional moment (M) as a function of non-dimensional time (t) with p,i 
varying from to 0.5 in increments of 0.1 for the degradation case. Values 



chosen were Ri/R 
m = 0.1. 
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Figure 2. (a) Non-dimensional moment (M) as a function of non- 
dimensional time (t) for various values of b\ varying from to 0.8 in increments 
of 0.2 for the degradation case. Here, Ri/R = 0.5, q — 1, ■0 = 1, D — 0.01, 
bo = n = 1, fii = 0.1, rii = 0.1. (b) Non-dimensional moment (M) as a 
function of non-dimensional time (t) for various values of ni starting at to 
0.8 in increments of 0.2 for the degradation case. Here, Ri/R = 0.5, q — 1, 
= 1, D = 0.01, b = n = 1, 6i = 0.1, /ii = 0.1. 
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Figure 3. Non-dimensional moment (M) as a function of non-dimensional 
time (t) for various values of the material parameters when the body is healing. 
(a) jli varying from to 0.5 in increments of 0.1 with bo = n Q = 1, bi = 0, 
n\ = (neo-Hookean model), (b) fii varying from to 0.5 in increments of 0.1 
with b = Uq = 1, &i = 0.1, rii = 0.1. (c) bi varying from to 0.8 in increments 
of 0.2 with bo = Uq = 1, Ji\ = 0.1, n\ = 0.1. (d) n\ varying from to 0.8 in 
increments of 0.2 with bo = no = 1, b\ = 0.1, Ji\ = 0.1. Other values chosen in 
(a), (b), (c), (d) were Ri/R = 0.5, q = l,$ = l,D = 0.01. 
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Figure 4. Non-dimensional angular displacement (ip) as a function of non- 
dimensional time (t) for various values of the material parameters when the 
body is degrading, (a) fli varying from to 0.5 in increments of 0.1 with 
bo = no = 1, bi = 0, rii = (neo-Hookean model), (b) p,± varying from to 
0.5 in increments of 0.1 with bo = no = 1, b\ = 0.1, n\ = 0.1. (c) b\ varying 
from to 0.8 in increments of 0.2 with bo = n = 1, fii = 0.1, n\ = 0.1. (d) n\ 
varying from to 0.4 in increments of 0.1 with 6 = ^o — 1> &i — 0.1, pi — 0.1. 
Other values chosen in (a), (b), (c), (d) were Ri/R = 0.5, q — 1, M — 0.5, 
5 = 0.01. 
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Figure 5. Non-dimensional angular displacement (ip) as a function of non- 
dimensional time (t) for various values of the material parameters when the 
body is healing, (a) Ji\ varying from to 0.5 in increments of 0.1 with bo = 
no = 1, b\ — 0, n\ — (neo-Hookean model), (b) fix varying from to 0.5 in 
increments of 0.1 with bo = no = 1, b\ = 0.1, n\ = 0.1. (c) b\ varying from 
to 0.8 in increments of 0.2 with b = n = 1, JJL\ = 0.1, n\ = 0.1. (d) n\ 
varying from to 0.4 in increments of 0.1 with 6 = n o — 1) &i — 0.1, pi = 0.1. 
Other values chosen in (a), (b), (c), (d) were Ri/R = 0.5, q — 1, M = 0.5, 
5 = 0.01. 
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Figure 6. (a) Non-dimensional moment (M) as a function of non- 
dimensional time (t) for various values of fti starting from to 0.5 in incre- 
ments of 0.1 for the degradation case with Ri/R = 0.75. (a) Non-dimensional 
moment (M) as a function of non-dimensional time (t) for various values of 
pi starting from to 0.5 in increments of 0.1 for the degradation case with 
Ri/R = 0.35. (c) Non-dimensional angular displacement ijp) as a function of 
non-dimensional time (t) for various values of fii starting from to 0.5 in incre- 
ments of 0.1 for the degradation case with Ri/R = 0.75. (d) Non-dimensional 
angular displacement (z/>) as a function of non-dimensional time (t) for various 
values of fix starting from to 0.5 in increments of 0.1 for the degradation case 
with Ri/R = 0.35. In all cases, q = 1, b = n = 1, &i = 0.0, n\ = 0.0. Initial 
and boundary conditions considered are given by 
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FIGURE 7. (a) Comparison of non-dimensional moment (M) as a function 
of non-dimensional time (t) for various constant diffusivity values for xjj — 1. 
(b) Comparison of non-dimensional angular displacement (ijj) as a function of 
non-dimensional time (t) for various constant diffusivity values for M — 1 . 
Also, Ri/R = 0.5, q = 1, b = n = 1, b\ = 0.1, rii = 0.1, /2i = 0.4. Initial 
and boundary conditions considered are given by 
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Figure 8. (a) Solution to the convection-diffusion equation (jH|) when the 
diffusivity is held constant, (b) Solution to (jUJ) when the diffusivity depends 
on the Almansi-Hamel strain, (c) Non-dimensional moment (M) as a function 
of non-dimensional time (t) with and without diffusivity depending on the 
Almansi-Hamel strain for ip = 1. (d) Non-dimensional angular displacement 
(■0) as a function of non-dimensional time (t) with and without diffusivity 
depending on the Almansi-Hamel strain for M — 1. Here, Ri/R = 0.5, q = 1, 
bo = no = 1, bi = 0.1, rii = 0.1, /2i = 0.4. For constant diffusivity case 
D = 0.01, and for diffusivity depending on the Almansi-strain, relation (|53|) is 
assumed with D = 0.01, Doo = 0.1, A = 1. Initial and boundary conditions 
considered are given by (145H47|) . 
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Figure 9. (a) Solution to the convection-diffusion equation ( jUJ) when the 
diffusivity is held constant, (b) Solution to (jUJ) when the diffusivity depends 
on the Almansi-Hamel strain, (c) Non-dimensional moment (M) as a function 
of non-dimensional time (t) with and without diffusivity depending on the 
Almansi-Hamel strain for ip = 1. (d) Non-dimensional displacement (ijj) as a 
function of non-dimensional time (t) with and without diffusivity depending on 
the Almansi-Hamel strain for M — 1. Here, Ri/R = 0.5, q = 1, b = n = 1, 
bi = 0.1, n\ = 0.1, fii = 0.4. For constant diffusivity case D = 0.01, and 
for diffusivity depending on the Almansi-strain, relation (|53|) is assumed with 
Do = 0.01, Ax, = 0.1, A = 1. Initial and boundary conditions considered are 
given by fl48 ti50|) . 
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Figure 10. Comparison of the steady state solutions (i = 40) to (|4"4"|) with 
constant diffusivity and when the diffusivity depends on the Almansi-Hamel 
strain, with initial and boundary conditions given by fj4*8H50|l . Ri/R = 0.5, 
q = 1, ip = 1, b = n = 1, b\ = 0.1, n\ = 0.1, pi = 0.4. For constant 
diffusivity case D = 0.01, and for diffusivity depending on the Almansi-strain, 
relation (j55|) is assumed with D = 0.01, D^ = 0.1, A = 1. 
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